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We investigate catalysis of electroweak baryogenesis by fermionic Higgs portal dark matter using 
a two Higgs doublet model augmented by vector-like fermions. The lightest neutral fermion mass 
eigenstate provides a viable dark matter candidate in the presence of a stabilizing symmetry Z 2 or 
gauged U(1 )d symmetry. Allowing for a non-vanishing CP-violating phase in the lowest-dimension 
Higgs portal dark matter interactions allows generation of the observed dark matter relic density 
while evading direct detection bounds. The same phase provides a source for electroweak baryogen¬ 
esis. We show that it is possible to obtain the observed abundances of visible and dark matter while 
satisfying present bounds from electric dipole moment (EDM) searches and direct detection experi¬ 
ments. Improving the present electron (neutron) EDM sensitivity by one (two) orders of magnitude 
would provide a conclusive test of this scenario. 


Introduction Precise cosmological observations have 
confirmed the existence of non-baryonic cold dark matter 
fl DMh 2 = 0.1187 db 0.0035 [1} [2]. How dark matter in¬ 
teracts with the Standard Model (SM) particles remains 
unknown. The discovery of the Higgs boson opens a new 
avenue for both probing the dark universe and modeling 
its dynamics. In particular, Higgs portal dark matter 0- 
[33] . which extends the SM with a SM gauge singlet that 
only interacts with the Higgs boson, provides one of the 
simplest realizations of this idea. 

It is possible that Higgs portal interactions may also 
provide an explanation for another unsolved problem 
in cosmology: the origin of the baryon asymmetry of 
the Universe (BAU). Combining data from Planck [I], 
WMAP [2j and large scale structure measurements, one 
has Yb = Pb/s = (8.59 =b 0.11) x 10 —11 , where pB is 
the baryon number density and s is the entropy density. 
Assuming that the Universe was matter-antimatter sym¬ 
metric at the end of the inflationary epoch, it is likely that 
interactions involving beyond the Standard Model fields 
generated the BAU during the subsequent cosmological 
evolution. To generate the observed BAU, these early 
universe interactions must satisfy the three Sakharov cri¬ 
teria [34]: (1) baryon number violation; (2) C and CP 
violation; (3) a departure from the thermal equilibrium 
(assuming CPT conservation). 

Fulfilling the criteria requires physics beyond the SM 
(BSM). Among the most attractive scenarios is elec¬ 
troweak baryogenesis p35H42] . as the relevant BSM in¬ 
teractions can be tested experimentally through direct 
searches for new particles at high energy colliders and 
probes of CP-violation (CPV) through searches for per¬ 
manent electric dipole moments (EDMs) [43j 04]. It is 
interesting to ask whether Higgs portal interactions re¬ 
sponsible for the observed dark matter relic density may 
also generate the BAU through electroweak baryogenesis 
(EWBG). 

In this paper, we explore this question with a simple 
Higgs portal realization involving fermionic dark matter. 


Denoting the latter by x, the lowest dimension Higgs 
portal interactions are non-renormalizable and have the 
form [23] 

, ( 1 ) 

where a , /? are effective couplings and A is the scale 
associated with the effective theory. Both interactions 
may yield an annihilation cross section consistent with 
the observed relic density under a thermal dark matter 
scenario. However, the first term with coefficient a is 
severely constrained by dark matter direct detection ex¬ 
periments XENON100 [45 j and LUX [46] . making it an 
unlikely candidate for dark sector physics. On the other 
hand, the sensitivity of direct detection searches to the in¬ 
teraction proportional to f3 is considerably weaker, due to 
the associated velocity suppression for a non-relativistic 
X- In what follows, we show that CPV Higgs portal in¬ 
teractions may lead to |/3| >> |cu| ~ 0, as needed for con¬ 
sistency with the observed relic density and direct detec¬ 
tion constraints, while also providing a sufficient source 
for the BAU through EWBG. 

Model We work in the Type-I two Higgs doublet 
model m~. wherein the second Higgs doublet has no 
Yukawa interactions with the SM fermions. By intro¬ 
ducing a pair vector-like fermions, 'ip l,r , transforming as 
(1, 2, 1/2) and Xl,r transforming as (1, 1, 0), the new 
Yukawa interaction and mass term can be written as 

YJlMffn + E (yM H iXR + y'iXlHUn ) + h.c. (2) 
/ i 

where / = -0, %. Clearly, fp^R may couple to the SM lep¬ 
tons through Yukawa interactions, precluding the possi¬ 
bility of a stable neutral fermion. To ensure stability, one 
may impose a Z 2 symmetry, in which new fermions are 
odd, while all the other fermions are even. Alternatively, 
one may also introduce a local U(1)d gauge symmetry, in 
which only ipL,R and Xl,r have non-vanishing but identi¬ 
cal charges. The new fermion mass matrix can be written 
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as 


( XL iPl) 


M x y'\Vi + y' 2 v 2 \ f Xr\ 
yivi + y 2 v 2 M.^ ) \tpR) 


(3) 


As observed in Ref. [35] the interactions © contain CP- 
violating phases that cannot be rotated way by field re¬ 
definition. 

The mass matrix in Eq. © can be diagonalized by per¬ 
forming separate flavor rotations on the left- and right- 
handed components of the new neutral fermions, viz 


Xl =Ml li £Li , 4 >l =Ul 2i (4) 


where £i ? 2 are the two mass eigenstates and where similar 
rotations are performed on the right-handed fields with 
matrix Ur. In the limit that these matrices are nearly 
diagonal, one has £i ~ X and £2 ~ V5 where -0 = VtL 
and x = xl + Xr- Assuming £1 « x is the lighter of the 
two states and does not decay on cosmologically relevant 
time scales (see below), one may obtain the effective dark 
matter-Higgs interactions by integrating out £2 and the 
charged component of the vector-like fermion doublet: 


X 


ViVi 


{cos 5iXX + i sin 5iXlsX} h\h^ -, 


(5) 


where = Arg (yiy'jM^), hj = (//,', //•' - v,), and the 
“H— •” denote interactions containing h\h 2 and (or) h\h\ 
or terms suppressed by x mixing. Comparing with 0 
we observe that a ~ cos Si and /3 ~ sin Si and A = M^. 
For sufficiently large | sin Si |, one may in principle satisfy 
the dark matter relic density and direct detection con¬ 
straints while providing a sufficiently strong CPV source 
for EWBG. 

Dark Matter For illustrative purposes, we assume the 
second Higgs doublet is superheavy, which is consistent 
with the current LHC measurements, so that the effec¬ 
tive interactions that dominate the annihilation of the 
DM are (x{l, 7s}x)^i^i- A more comprehensive study 
of the dark matter phenomenology in the 2HDM will be 
given a forthcoming paper [49]. For the annihilation in 
the early universe responsible for setting the thermal relic 
density, dark matter is moderately relativistic. For direct 
detection, on the other hand, the dark matter particles 
are non-relativistic, with typical velocity v ~ 10 -3 in the 
galactic halo m • As a result, the parity-violating ef¬ 
fective coupling contributes to dark matter annihilation 
but gives a velocity suppressed contribution to the dark 
matter-nucleus scattering cross section. It is, thus, possi¬ 
ble to choose values of sin Si/M/ that allow for sufficiently 
efficient annihilation as needed to obtain the observed 
relic density while evading the stringent spin-independent 
direct detection bounds. 

We illustrate this feature in Fig. [l] choosing y\ = ye 16 
with S being the only CP phase, y[ = y 2 =0.5 , y f 2 = 0, 
tan f3 = 15, M x = 300 GeV and M ^ = 1500 GeV. In this 


case the heavy-light (x-VO mixing is 0(1O -2 ) and can be 
safely neglected. The resulting region between the red 
dashed lines gives a dark matter relic density consistent 
with Planck and WMAP measurements at 95% C.L. . 
The grey hatched region leads to an over-abundance of 
DM, while the area to the upper right of the red-dashed 
curves corresponds to under-saturating the relic density. 
Regions below the black solid line are excluded by the 
LUX result [46] . 

Electric dipole moments The CP-odd Yukawa cou¬ 
plings given in © generate an elementary EDM via two- 
loop Barr-Zee diagrams m- For illustrative purposes, 
we will work in the limit that the masses of the remain¬ 
ing neutral scalars and charged scalars are sufficiently 
heavy that the dominant contributions arise from the ex¬ 
change of a W + W~ pair. CP-violation enters the con¬ 
tribution through the relative phase of left- and right- 
handed charged currents. A direct calculation gives the 
following result for the electron EDM 

d e « 4 2;)s (v Im (^L2^R2i)^U^/vvvi/(n,r2) (6) 

m w 

where d^ ~ 2.5 x 10 -27 e • cm, sw = sin 0w with Qw 
being the weak mixing angle, rq = (m^+ /mw) 2 and r 2 = 
( mi/mw ) 2 , and the loop function fww{%,y) is given in 
Ref. [52]. An analogous expression applies for the quark 
EDMs that give rise to the EDM of the neutron. 

Present constraints obtained by the ACME collabo¬ 
ration [53] yield the most stringent limit on the elec¬ 
tron EDM: d e = (-2.1 d= 4.5) x 1(T 29 e cm . This 
bound is nearly two orders of magnitude more stringent 
than the current neutron EDM limit [54]. However, given 
the order-of-magnitude larger value for the light quark 
Yukawa couplings, the electron and neutron EDM con¬ 
straints on this scenario differ by only by roughly a factor 
of ten. Nonetheless, the bounds from d e are presently 
too weak to impact the parameter space shown in Fig. [l] 
largely because only the W + W~ graph contributes to 
the EDM (the situation is analogous to the two-loop 
chargino-neutralino EDM contributions in “split super- 
symmetry” [52l [55l [56] 5. Planned future improvements 
in EDM search sensitivities would change this situation. 
To that end, we show in Fig. [l]the parameter space that 
would be probed by a possible future d e search that is 
an order of magnitude more sensitive than the present 
ACME bound (see [48] for a detailed analysis). Regions 
to the right of blue dotted line would be excluded should 
a null result be obtained with this sensitivity. The inputs 
are the same as these in making the dark matter curves. 
Baryogenesis The three Sakharov conditions are real¬ 
ized in the following way. First, a strongly first order 
electroweak phase transition, which provides a departure 
from equilibrium [57H59] , is induced by the two Higgs 
doublet potential. The phase transition occurs as bubbles 
of broken electroweak symmetry, which nucleate and ex¬ 
pand in a background symmetric phase, fill the Universe. 
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Second, CP violation arises from the Yukawa couplings 
of the new fermions. A non-zero CP-asymmetric charge 
density is produced via the CP-violating interactions at 
the walls of the expanding bubbles, where the Higgs vac¬ 
uum expectation value (VEV) is space-time dependent. 
This charge density diffuses ahead the expanding bub¬ 
bles and is transformed into charge densities of other 
species via inelastic interactions in the plasma, result¬ 
ing in a net density of left-handed fermions Ur. Third, 
the weak sphaleron processes, which are unsuppressed in 
the symmetric phase, violate the baryon number. The 
presence of nonzero tir biases the sphaleron processes, 
leading to the production of the baryon asymmetry m- 
The expanding bubbles capture the asymmetry, which is 
preserved assuming the first order EWPT is sufficiently 
strong so as to quench the sphaleron transitions inside 
the bubbles. 

We compute the CP-asymmetry and resulting ur using 
the closed-time path formalism and the “VEV-insert ion” 
approximation [61, 62,. Under these assumptions, one 
obtains a coupled set of diffusion equations 

« = -£ r A(/*-Mi--)+ ^ CPV (7) 

Aj 


where j£ is the number density current for species of 
particle k ; /i&j are chemical potentials for species &, j 
etc ; the Fa denote rate for particle number changing re¬ 
actions; and S£ PV is the CP-violating source term asso¬ 
ciated with the CP-violating phases and the space-time 
varying vacuum expectation values of the neutral Higgs 
scalars. 

In the present instance, the CP-violating source enters 
the diffusion equations for the new neutral fermions x 
and the neutral component of \/j. Details of the compu¬ 
tation for the model under consideration here are similar 
to those given in see Ref. [48], so we do not repeat all of 
them here. As in obtaining the other curves in Fig. [l] we 
assume y\ contains the only CP phase and y' 2 = 0. The 
source S^ PV = —5^ PV depends crucially on the damping 
rate of the new fermions, corresponding to the imaginary 
parts of the poles in the finite temperature fermion prop¬ 
agators [63HH5] : 


7^ 


3g 2 + g' 2 + 4a 2 ^ 
32 t r 


7x 


a2 9 2 N T 

Stt 


(8) 


where gjy is the gauge coupling of the U(1)d and a is the 
relevant hyper-charge of ^ and y. We have neglected the 
contribution of the Higgs mediated processes, which is 
subdominant compared with the contribution from scat¬ 
tering with gauge bosons [64]. 

The expression in Eq. © is general, and one must 
apply it to the specific representations of interest here. 
For the i/j, one must include the contributions from the 
SU(2U(l)y, and U(1 )e> gauge bosons when the lat¬ 
ter are included in the model to provide for dark matter 


stability; for the x one need only include the new U(1)d 
contributions (Higgs mediated interactions contribute to 
both). Notice that a definitive quantitative treatment of 
the CPV fermion sources remains an open problem. The 
VEV insertion approximation used in our calculation is 
likely to overly estimate the CPV source by at least a 
factor of a few. The results given here, thus, provide 
a conservative basis for evaluating the dark matter and 
EDM restrictions on the parameter space viable for the 
EWBG. We refer to [35] and references therein for a de¬ 
tailed discussion of the theoretical issues associated with 
the calculation of the CPV source term. 

Solution of the coupled set of diffusion equations ([7]) 
can be simplified through several observations. First, 
since all light quarks have nearly equivalent diffusion 
constants and are mainly produced by strong sphaleron 
processes, baryon number conservation on time scales 
shorter than the inverse of the electroweak sphaleron rate 
approximately implies the relations qiL = q 2 L = —2ur = 
—2 cIr = —2 sr = —2 cr = —2 bR —2b = 2(Q + T). The 
resulting set of equations for new fermions are already 
given in Ref. [48]. 

Second, one may obtain a set of equations solely for 
the number densities under the diffusion approximation: 
J = — D\/n , where D is the diffusion constant for a given 
species whose number density three-current is J. The 
diffusion constants for quarks, leptons, and Higgs scalars 
have been computed previously [66]. Here we compute 
those for the new fermions, focusing elastic scattering, t- 
channel vector boson and Higgs scalar exchange diagrams 
that are expected to dominate the scattering process and 
play the major role in limiting the diffusion of particles. 
The diffusion constants mediated by the gauge boson and 
Higgs bosons can then be written as 


Dp) 

D t 


60 
7lT 


a 2 D a 2 T In 
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7n 167t 2 
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( 10 ) 


The inverse diffusion constants for x and V’ are then: 
D- 1 * D~ b ) + E,; D h ) and D^ 1 * T /100 + D B } + 
E,: D h) > where M§, is the one-loop Debye mass of 
the new gauge boson, M\, = ^ a 2 D T 2 [66] |67], where 
cxr> = Mjj. are the thermal mass of the “i”th 

Higgs boson. 

We solve diffusion equations analytically by assuming 
that strong sphaleron transitions as well as the Yukawa 
interactions for the top quarks and the new fermions 
are all in chemical equilibrium. These assumptions lead 
to a series of conditions relating the chemical potentials 
and, therefore, the number densities of various species. 
Chemical equilibrium for the new fermion Yukawa inter¬ 
actions implies that n x /k x — nn/kH —n^/k^ = 0, where 
rii and ki {i = y, ^ H) are the number density and 
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the statistical factor for particle “i” respectively. Fur¬ 
thermore, in the static limit ( where v w —)> 0), we have 
D x n x = where we have assumed the boundary 

conditions n x ( oc) = n' x ( oc) = n^(oo) = n^(oc) = 0. 
Therefore, all the charge densities can be expressed as 
the functions of the Higgs density, and all that remains 
is to solve for the Higgs density. We refer to Refs. [68j 69 
for the semi-analytical approximate solution in this case. 

The left-handed charge density that biases weak 
sphaleron transitions, is given by the sum of all charge 
densities of left-handed quarks and leptons of all gen¬ 
erations, til = EM + &i)- The left-handed density 
is converted into a baryon number density n B through 
the weak sphaleron process. The following formula de¬ 
scribes baryon generation and washout ahead of the bub¬ 
ble wall [70] 

r^s f 7 / \ i 

n B = -n F - — / dz n L (z)e 4 (11) 

^V w J — oo 

where v w is the bubble wall velocity, T ws is the weak 
sphaleron rate and 2: is the spatial coordinate perpendic¬ 
ular to the wall in the frame where the wall is at rest. 
Negative values of z correspond to the symmetric elec- 
troweak phase, positive values to the broken phase. 

For explicit numerical illustration, we consider the case 
wherein stability of the lightest neutral new fermion is 
guaranteed by the U{1) B symmetry. In addition, we take 
the Higgs vev profiles in Ref. [M, [72] and follow Ref. m 
E2] in our calculation of the CP-conserving relaxation 
rates for quarks and Higgs. We assume &Q3 = 2 k F = 
2 k B = 6, k B = 4, k^ = 2 k x = 2 and g n = 0.35 for 
numerical calculation. 

In Fig. [TJ we show contours of constant Y B in the y — S 
plane. Regions between the green dot-dashed lines are 
consistent with the observed Y B at 95% C.L.. As noted 
earlier, the current ACME result places no constraint on 
the CP-violating phase, so one has adequate parameter 
space for successful baryogenesis. However a prospective 
future search for the electron EDM could probe most of 
the parameter space. A similar statement would apply to 
a future neutron EDM search, assuming an improvement 
in sensitivity of two orders of magnitude. Since S domi¬ 
nates the strength of the parity-violating DM-Higgs cou¬ 
pling, the present exclusion limit given by the LUX puts a 
strong lower bound on the phase. Nevertheless, it is pos¬ 
sible to find parameter space where both sizable baryon 
asymmetry and a successful dark matter candidate can 
be derived. We note that the dark matter/EWBG-viable 
parameter space varies significantly as one changes inputs 
of M x and M^. A much smaller mass splitting between 
M x and M ^ would induce an enhanced source term, via 
the resonant enhancement [61], for the EWBG, but in¬ 
troduces a large mixing between two neutral fermions, 
which may degrade the viability of y as a dark matter 
candidate. A much larger mass splitting would lead to 
an inadequate baryon asymmetry on the other hand. 



y 


FIG. 1: Contours of constant baryon asymmetry Yb in units 
of 10 -10 in the y — 5 plane (see text for parameter definitions). 
Dashed green band indicates a value of Y B consistent with 
the observed asymmetry, while the red dashed band indicates 
parameter region consistent the observed dark matter relic 
density. Grey hatched region leads to an over abundance of 
dark matter. Dark matter direct detection constraints from 
the LUX experiment m allow the region above the solid 
black line. A future search for the electron EDM with an order 
of magnitude greater sensitivity than the ACME result [53] 
would probe the parameter space above the blue dashed line. 

Conclusion The possibility of correlating the abundance 
of visible and dark matter is a topic of considerable in¬ 
terest, inspiring new paradigms such as asymmetric dark 
matter [73] (for recent reviews, see Refs. mi usd or 
WIMP baryogenesis m- In this study, we have investi¬ 
gated the possibility that CP-violating, fermionic Higgs 
portal interactions may both govern the abundance of 
the dark matter relic density and drive the production 
of the baryon asymmetry through EWBG while evad¬ 
ing the present EDM and dark matter direct detection 
bounds. Working with a simple two Higgs doublet model 
augmented by the presence of vector-like fermions, we il¬ 
lustrated the existence of parameter space regions con¬ 
sistent with the observed abundance of visible and dark 
matter and present experimental constraints. An exten¬ 
sive study of the dark matter /Y B -viable parameter space, 
along with the relevant phenomenological implications, 
will appear in aforthcoming publication. 
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